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1.  INTRODUCTION1 


The  objective  of  this  study  is  to  develop  the  necessary  techniques  and 
computer  codes  to  theoretically  determine  the  effects  of  rocket  exhaust  plumes 

on  radio  frequency  (RF)  signal  propagation  between  the  launch  missile  and  a 
ground  or  airborne  station.  Primary  emphasis  will  be  placed  on  the  altitude 
regime  for  the  third  stage  solid  engine  of  the  MX  missile.  The  focus  of  the 
studies  conducted  thus  far  were  to  review  various  aspects  of  the  state  of  the 
art  in  RF  interference  from  rocket  plumes,  to  study  and  determine  c he  dominant 
factors  in  affecting  the  plume  plasma  and  in  influencing  the  RF  plume  inter¬ 
action,  and  to  define  the  state-of-the-art  models  required  for  its  solution. 
The  results  and  findings  are  presented  in  this  report. 


Continuous  communication  with  a  missile  during  launch  is  essential  for 
the  successful  and  safe  completion  of  the  Launch  operation.  Critical  opera¬ 
tions,  such  as  target  tracking,  telemetry  transmission,  and  command-destruct , 
all  depend  upon  adequate  propagation  of  electromagnetic  radiation  in  the  radio 
wave  spectrum  from  or  to  the  launch  vehicle.  Unfortunately,  it  often  occurs 
that  the  rocket  exhaust  may  be  in  a  high  level  of  ionization  which  gives  rise 
to  free  electrons  in  the  exhaust  plume.  These  free  electrons  can,  in  many 
cases,  interfere  with  the  propagation  of  an  electromagnetic  signal  when  they 
are  present  in  sufficiently  large  amounts  and  when  the  line  of  sight  between 
the  missile  and  the  ground  antennas  traverses  the  missile  plume. 


The  amount  of  electrons  in  the  plume  is  intimately  related  to  the  alkali 
metal  impurities  in  the  propellant,  for  the  rocket  engine  temperature  is  gener¬ 
ally  not  high  enough  to  cause  ionization  of  the  other  combustion  products. 
However,  because  the  alkali  metals  (notably  potassium  and  sodium)  have  low 
ionization  potentials,  they  become  highly  ionized  with  the  result  that  consid¬ 
erable  amounts  of  electrons  are  generated  and  deposited  in  the  exhaust  plume. 
In  general,  solid  rocket  propellants  contain  much  larger  amounts  of  sodium  and 
potassium  contaminants  than  liquid  propellants.  This  is  mainly  due  to  the 
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chemical  process  of  manufacturing  ammonium  perchlorate,  an  essential  ingredient 
in  the  solid  propellant,  which  carries  over  some  original  raw  materials  con¬ 
taining  alkali  metal  compounds.  While  this  study  is  primarily  aimed  at  exnaust 
»-%  i it>  cj  f  '**  o,Ti  5  0^  ^  d  **  p  C  *  nt  ^ '  n  »  s  j  ti  !*i  ^  .d  n  cj  1 v  i.  *5  .1  ”  ?•  L s  js  t  o  £  ^  o  *2  r  r L  *.  ■? 1 1  t.  * 
rockec  engines  ns  well*  The  principal  differences  between  t’ne  two  types  oc 
rocket  plumes  lie  in  the  presence  of  large  amouncs  of  particulates  in  the  solid 
rocket  plume  and  the  different  chemical  composition.  Therefore,  the  ability  to 
calculate  solid  rocket  plumes  automatically  includes  liquid  rocket  plume  by 
deleting  the  particular  phase  in  the  plume  flow  field. 


RF  interferences  from  rocket  plumes  have  been  observed  in  numerous 
launch  operations  in  the  past.  References  1  and  2  provide  good  summaries  of 

the  existing  data.  More  detailed  description  of  field  data  from  various  launch 
sites  may  be  found  in  Refs.  3  through  9.  A  large  amount  of  S-band  telemetry 
attenuation  data  are  also  available  from  numerous  Minuteman  III  launches  from 


^Victor,  A.C . j  "Plume  Signal  Interference,  Part  1:  Radar  Attenuation,"  Report 
NWC-TP-5319,  Part  1,  Naval  Weapons  Center,  China  Lake,  CA,  June  1985. 

2Victor,  A.C.,  "Plume  Signal  Interference,  Part  1:  Plume-Induced  Noise  (U)," 
Report  NWC-TP-5319,  Part  2,  Naval  Weapons  Center,  China  Lake,  CA,  June  1975 
(Confidential) . 

2Poehler,  H.A. ,  "Results  of  Special  Flame  Measurements  of  Titan  III-C  Tests 
3656/6025  and  6020/6546,"  Report  ETV-TRM-66-30 ,  Pan  American  World  Airways, 
Guided  Missile  Range  Division,  1966. 

^Poehler,  H.A.,  "Project  See-Thru,  Flame  Attenuation  Interference  Measure¬ 
ments,  Titan  III-C  Test  8275/2250,  Preliminary  Report,"  Report  No.  ETR-TR- 
68-6,  Pan  American  World  Airways,  1967. 

^Poehler,  H.A. ,  "Project  See-Thru,  Flame  Interference  Measurements,  Titan 
III-C  Launch  Test  8275/2250,  Final  Report,"  Report  No.  ETR-TR-68-3,  Pan 
American  World  Airways,  1968. 

6Poehler,  H.A. ,  "Rocket  Exhaust  Signal  Attenuation  and  Degradation,  Final 
Report,  Vols.  I,  II,  III,"  Report  No.  ETR-TR-69-4,  Pan  American  World  Airways, 
Aerospace  Service  Division,  1969- 

7Eiy,  O.P.,  "Rocket  Exhaust  Effects  on  Radio  Frequency  Transmission," 

J.  Spacecraft,  PP*  310-314,  1966. 

sVicente,  F.A.,  E.C.  Taylor,  and  R.W.  Phelps,  "Analysis  of  Flame  Effects  on 
Measuring  Electromagnetic  Propagation  Data,"  J.  Spacecraft,  4,  pp.  1069-1075, 
1967. 

^Golden  K.E.,  E.L.  Taylor,  and  F.A.  Vicente,  "Diffraction  by  Rocket  Exhaust," 
Communication  in  IEEE  Transaction  on  Antenna  and  Propagation,  AP16 , 
pp.  614-616,  1968. 
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Western  Test  Range  (WTR).  They  were  recorded  at  antenna  locations  over  some 
distances  apart  and,  therefore,  should  provide  indications  or  the  aspect  angle 
variations  of  plume  interference.  They  are  available  as  Post-Launch  Analysis 
of  Telemetry  Systems  (PLATS')  reports  issued  by  Western  Space  and  Missile 
Center  (WSMC). 10,11 

These  data  have  shown  both  amplitude  and  phase  interference  of  the  RF 

signal  due  to  the  presence  of  the  rocket  exhaust  plume.  In  many  instances, 

both  amplitude  and  phase  noise  were  also  detected.  It  appears  that  the  noise 

may  be  due  to  turbulent  fluctuations  in  the  plume  plasma.  Spectral  analysis  of 

12 

the  amplitude  noise  has  been  made  and  indicated  a  peak  noise  at  about  1  kHz. 

,'lthough  it  is  necessary  to  account  for  amplitude  and  phase  perturba¬ 
tions  in  the  RF  signal  in  the  mean  values  as  well  as  noises  to  completely 
define  rocket  plume  effects  on  RF  signal  propagation,  the  initial  efforts  of 
this  will  be  confined  to  determining  the  effects  on  the  mean  signal 

arplitude  a^d  phase.  Only  the  mean  plume  flow  field  properties  will  be  calcu- 
Jared,  ana  the  RF-plasma  interaction  analysis  will  be  restricted  to  employing 
th?  mean  plasma  properties  and  determining  the  mean  amplitude  and  phase 
changes.  This  is  necessary  in  order  to  confine  the  initial  efforts  to  a 
tractable  level.  The  inclusion  of  amplitude  and  phase  noise  study  will  neces¬ 
sitate  treating  turbulent  fluctuations  in  the  rocket  plume  flow  field  and  may 
entail  studies  beyond  the  present  state  of  the  art.  It  appears  that  at  the 
altitudes  of  interest  associated  with  the  third  stage  burn  of  the  MX  missile, 
the  turbulence  in  the  rocket  plume  is  mainly  originated  from  the  combustion 


10"Post  Launch  Analysis  Telemetry  Systems  Preview — Operation  3162,"  Report  No. 
PA300-T-79-31P ,  Systems  Performance  Analysis  Department,  Western  Test  Range 
Division,  Federal  Electric  Corporation,  26  July  1979. 

^"Post  Launch  Analysis  of  Telemetry  Systems — Operation  5544,  Final,"  Report 
No.  A.S300-T-81-21F,  Systems  Performance  Analysis  Department,  Western  Test 
Range  Division,  Federal  Electric  Corporation,  1  April  1981. 

^McKelvey,  G.R.,  "Exhaust  Plume  Test  Results — Report  Number  2,"  Technical  Note 
PA300-N-79-23 ,  Systems  Performance  Analysis  Department,  Western  Test  Range 
Division,  Federal  Electric  Corporation,  27  September  1979. 
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chamber  of  the  rocket  engine  and  retained  by  the  core  flow  of  the  plume.  This 
is  in  contrast  to  plume  turbulence  at  low  altitudes,  such  as  during  the  first 
stage  burn,  where  t'na  dominant  plume  turbulence  Ls  generated  by  che  free  shear 


layer  between  the  plume  gas  and  the  ambient  air.  While  a  li 

i  3 

data  of  turbulent  fluctuations  in  the  plume  core  does  exist,' 


ted  amount  c : 
it  pertains 


to  those  from  a  liquid  rocket  engine  of  small  thrust.  Applications  to  solid 


rocket  engines  and  higher  thrust  levels  require  further  study.  The  inclusion 
of  plume  turbulence  is  beyond  the  scope  of  the  present  initial  effort. 


The  results  and  conclusions  of  this  report  consist  of  various  studies 

and  analyses  to  define  the  necessary  theoretical  models  and  computer  codes 
that  determine  the  RF  signal  amplitude  and  phase  interference  due  to  rocket 
p lumes . 


^Sutton,  G.A. ,  "Boost  Measurements  and  Analysis  Program  Final  Report,  Vol.  7: 
Plasma  Properties  Derived  From  Full  Scale  Rocket  Motor  Measurements,"  Ford 
Aerospace  and  Communications  Corporation  Report  No.  U-6511,  prepared  by 
Concord  Sciences  Corporation,  Concord,  MS,  15  March  1979. 

^August  G.  and  D.  Tremain,  "Boost  Measurements  and  Analysis  Program  Final 
Report,  Vol.  9:  Rocket  Plume  Measurements  and  Data  Evaluation,"  Ford  Aero¬ 
space  and  Communications  Corporation  Report  No.  U-6511,  prepared  by  Concord 
Sciences  Corporation,  Concord,  MS,  15  March  1979. 


12 


2.  METHODOLOGY  OF  RF-PLUME  PLASMA  INTERACTION  ANALYSIS 


The  analyses  necessary  for  determining  the  plume  effects  on  AF  signal 
propagation  are  complex  and  lengthy.  It  is  useful  to  separate  the  analyses 
into  a  number  of  sequential  calculations  for  various  chemical  and  physical 
phenomena : 

o  Rocket  Engine  Characterization  and  Exit  Plane  Flow  Properties 

o  Rocket  Plume  Flow  Field 

o  RF-Plutne  Plasma  Interaction 

These  calculations  are  briefly  described  below. 

2.1  ROCKET  ENGINE  CHARACTERIZATION  AMD  EXIT  PLUME  FLOW  PROPERTIES 

The  first  step  in  the  analysis  is  to  determine  the  properties  of  the 
rocket  engine  exhaust  at  the  nozzle  exit-  The  essential  data  to  carry  out 
this  calculation  are  the  chemical  composition  of  the  rocket  propellant 
including  the  amount  of  alkali  metal  contaminants,  the  rocket  chamber 
pressure,  and  the  detailed  geometrical  configuration  of  the  engine  chamber, 
nozzle  throat,  and  nozzle  contour. 

The  availability  of  the  alkali  metal  contaminant  levels  in  the  pro,,?.l- 
lant  is  essential  to  the  determination  of  the  electron  density  in  the  rocket 
plume.  These  contaminants  are  introduced  in  the  propellant  mainly  via  ammo¬ 
nium  perchlorate,  an  essential  ingredient  in  the  common  solid  propellants.  A 
rule  of  thumb  to  estimate  their  content  is  to  assume  200  parts  per  million 
(ppm)  of  potassium  and  125  ppm  of  sodium  by  weight  to  be  present  in  the  ammo¬ 
nium  perchlorate*  and  then  to  compute  their  weight  fraction  in  the  composite 


*Private  communication  with  E.M.  Landsbaum,  The  Aerospace  Corporation. 
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propellant.  This  is  only  a  very  crude  estimate,  however,  and  it  often  produces 
unreliable  results.  It  is  imperative  that  the  propellant  be  chemically  ana¬ 
lysed  for  alkali  imourity  content  to  define  the  weight  fraction  accurately. 


i.ne  characterisation  of  tne  rocaet  e xr, a u s c  tiow  properties  begins  vita  a 
calculation  of  the  rocket  engine  chamber  flow  properties  and  composition  by 
assuming  chemical  equilibrium  to  prevail  at  the  specified  chamber  pressure  at  a 
constant  enthalpy.  3ecause  the  flow  velocity  is  low  and  the  chamber  tempera¬ 
ture  is  high,  the  combustion  time  is  very  short  in  comparison  with  the  time 
required  for  the  combustion  products  co  flow  through  the  entire  chamber.  Under 
this  condition,  chemical  equilibrium  is  attained. 


As  the  combustion  products  flow  into  the  entrance  region  of  the  nozzle 
throat,  their  velocities  increase  rapidly  and  deviation  from  chemical  equili¬ 
brium  starts.  Therefore,  the  calculation  must  consider  finite  chemical  reac¬ 
tion  rates  as  the  engine  exhaust  flows  into  the-  thnat  and  is  expanded  in  the 
rocket  nozzle.  Since  a. large  amount  of  particulates  is  present  in  the  combus¬ 
tion  products  of  solid  propellants,  it  is  often  necessary  to  treat  two-phase 
flow  in  the  throat  and  nozzle  calculations. 


While  a  number  of  computer  codes  exist  for  the  calculations  described 
in  this  section,  the  Solid  Performance  Program  (SPP)  code  developed  under  the 
auspices  of  the  Air  Force  Rocket  Propulsion  Laboratory  combined  them  into  one 
code  which  has  been  adapted  for  the  present  study. 

The  gas  velocity  is  subsonic  in  the  chamber  and  becomes  supersonic 
downstream  of  the  nozzle  throat.  Much  attention  has  been  devoted  to  the  method 
of  computing  the  transonic  flow  as  the  two-phase  rocket  exhaust  passes  from  the 
chamber  to  the  nozzle.  The  approximate  transonic  analysis  of  the  SPP  code 
avoids  treating  elliptic  partial  differential  equations  for  the  subsonic  flow 
by  solving  the  inlet  region  of  the  nozzle  as  a  one-dimensional  flow  problem. 
As  the  exhaust  passes  through  the  throat,  an  uncoupled  particle-gas  flow  with 
constant  fractional  thermal  a".d  velocity  lag  was  assumed.  Solutions  from  the 
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SPP  code  will  be  compared  with  the  cylindrical  three-dimensional  two-phase  flow 
(CY3D2P)  code“^  developed  at  The  Aerospace  Corporation  to  check  the  accuracy  of 
the  former  code.  The  CY3D2?  cole  treats  the  fully  coupled,  unsteady  governing 
equations  for  both  ch_  gas  and  the  particulate  phases  so  as  to  :  i  tain  the 
hyperbolic  form  of  the  governing  equations.  It  provides  an  exact  solution  to 
the  two-phase  flow  problem;  however,  it  does  so  ac  the  expense  of  additional 
computing  time  and  storage. 


Typical  results  from  the  SPP  code  for  the  third  stage  of  the  MX  missile 
as  well  as  its  comparison  with  Che  CY3D2P  code  will  be  presented  Later  in  this 

report.  It  is  anticipated  that  when  the  SPP  code  is  thus  validated,  it  will 
be  adopted  for  all  subsequent  calculations. 


ROCKET  PLUME  FLOW  FIELD 


The  next  step  in  the  analysis  is  to  take  the  rocket  engine  exhaust  flow 

properties  at  the  nozzle  exit  plane  to  continue  the  calculation  outside  the 

nozzle.  When  only  gas  phase  is  considered,  a  number  of  simplified  approxima- 

16  —  18 

tions  may  be  used  for  the  far  field  of  the  rocket  plume.  The  approxima¬ 

tions  assume  chat  the  gas  has  been  sufficiently  expanded  outside  the  nozzle  to 
approach  a  source-like  flow,  and  the  detailed  mechanisms  of  the  gas  expansion 
in  the  vicinity  of  the  nozzle  are  of  minor  importance.  Their  ability  to 
describe  the  plume  gas  density  distribution  was  investigated.  Also,  in  the 
low  altitude  regime  of  the  MX  third  stage  burn  for  the  B0A3  trajectory,  inter¬ 
action  between  the  plume  gas  and  the  ambient  atmosphere  becomes  important  in 


^Chang,  1-3.,  "Three-Dimensional,  Two-Phase  Supersonic  Nozzle  Flows,"  AIAA 
Paper  No.  81-1219,  AIAA  14th  Fluid  and  Plasma  Dynamics  Conference,  Palo  Alto, 
CA,  23  June  1981. 

^HiLl,  A.F.J.  and  J.S.  Draper,  "Analytical  Approximation  for  the  Flow  from  a 
Nozzle  into  a  Vacuum,"  J.  Spacecraft,  pp.  1552-1554,  1966. 

^Brook  J.W.,  "Far  Field  Approximation  for  a  Nozzle  Exhausting  into  a  Vacuum," 
J.  Spacecraft,  6,  pp.  626-628,  1969. 

^Boynton,  F.P.,  "Exhaust  Plumes  from  Nozzles  with  Wall  Boundary  Layers," 

J.  Spacecraft,  5.,  pp.  1143-1147,  1968. 
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affecting  the  plume  plasma  properties.  The  interaction  generates  the  bow  shock 

and  the  barrel  shock  in  the  plume.  Its  effect  on  the  plume  shape  was  consid- 

19  . 

ered.  The  method  of  characteristics  (MOC)  solution  has  been  applied  to  tne 
inviscid  flow  of  plume  in  the  past  with  some  success  and  was  used  in  this  -icu-iy 
to  compute  the  pLume  flow  field  and  pLasma  properties.  It  has  t h ->  capability 
of  treating  the  barrel  shock  formation  while  assuming  modified  Newtonian  pres¬ 
sure  at  the  plume  boundary. 


When  the  presence  of  particulates  in  the  plume  needs  to  be  considered, 

the  computation  becomes  more  complex  and  lengthy.  Two  computer  codes  are 

available  for  calculating  two-phase  plume  flow  fields  by  the  method  of  charac- 

20  21 

teristics,  namely,  the  Reacting  Multi-Phase  Code  (RAMP)*-  ’  and  the  Two- 

22  23 

Dimensional  Two-Phase  Plume  Analysis  Program  (TD2P).  ’  The  latter  code  was 

used  in  chis  work  to  study  the  effects  of  particulates  on  the  plume  gas  density 
distribution.  The  effects  of  the  presence  of  the  particulates  will  be  dis¬ 
cussed  in  Section  3.3.2. 

More  ionization  in  addition  to  that  occurring  in  the  nozzle  was  found  by 

24  ... 

the  Martin-Marietta  Corporation  when  the  region  just  outside  the  nozzle  was 


^Sraith  S.D.  and  A.W.  Tatliff,  "Rocket  Exhaust  Plume  Computer  Program  Improve¬ 
ment,  Volume  IV:  Final  Report,"  Report  No.  LMSC/HREC  D162220-IV-A,  Lockheed 
Missiles  and  Space  Company,  Huntsville,  AL,  June  1971. 

2^Penny,  M.M. ,  et  al.,  "Supersonic  Flow  of  Chemically  Reacting  Gas-Particle 
Mixtures,  Volume  I:  A  Theoretical  Analysis  and  Development  of  the  Numerical 
Solution,"  Report  No.  LMSC/HREC  TR  D496555-I,  Lockheed  Missiles  and  Space 
Company,  Huntsville,  AL,  January  1976. 

2*Penny,  M.M.,  et  al.,  "Supersonic  Flow  of  Chemically  Reacting  Gas-Particle 
Mixtures,  Volume  II:  RAMP — A  Computer  Code  for  Analysis  of  Chemically  React¬ 
ing  Gas-Particle  Flows,"  Report  No.  LMSC/HREC  TR  D496555-I..,  Lockheed 
Missiles  and  Space  Company,  Huntsville,  AL,  January  1976. 

22Kliegel,  J.R.  and  G.R.  Nickerson,  "Axisymmetric  Two-Phase  Perfect  Gas  Perfor¬ 
mance  Program,  Volume  I:  Engineering  and  Program  Description,"  Report  No. 
02784-6006-R000,  TRW,  3  April  1967. 

2J"User's  Manual  for  the  Two-Dimensional  Two-Phase  Plume  Analysis  Computer 
Program,"  Report  No.  SN-141,  Dynamic  Sciences,  Inc.,  December  1968. 

2^Roraine,  G.L.,  et  al.,  "RF  Signal  Degradation  from  the  MX  Stage  III  Plume," 
Report  No.  SE4-116037,  Martin-Marietta  Corporation,  Denver,  CO,  May  1981. 
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considered.  Presumably,  their  calculation  was  done  in  the  region  near  the 
thrust  axis  where  the  expansion  tan  from  the  nozzle  tip  does  not  reach  until 
further  downstream.  sic’ever ,  near  the  no.::'.e  lip,  the  expansion  is  much  more 
r  »pid  and  should  cause  the  chemical  compos  II  inn  to  freeze  ac  the  exit  val..:. 
Thus,  downstream  of  the  exit  plane,  variations  of  the  chemical  composition  may 
exist.  This  aspect  of  the  additional  ionization  has  not  been  fully  investi¬ 
gated  and  would  requi.e  additional  ana’ysis.  At  present,  we  assume  the  elec¬ 
tron  mole  fraction  in  the  whole  pluiTj  to  be  at  the  nozzle  exit  value.  The 
electron  density  anywhere  in  the  plume  can  then  be  computed  from  the  gas 
density  distribution.  Together  with  the  plume  temperature  distribution,  if 
enables  the  collision  requency  to  be  calculated.  These  are  the  plasma  prop¬ 
erties  required  for  de  ermining  the  interference  of  the  electromagnetic  waves 
by  the  rocket  plume. 

2.3  RF-PLl/ME  PLASMA  INTERACTION 


Analysis  of  the  . interac t ion  between  the  RF  signals  and  the  plasma  is 

based  on  the  concept  of  an  effective  dielectric  constant,  e^,  which  depends  on 

the  ratios  of  plasma  frequency  and  collision  frequency  to  the  signal  fre- 
25 

quency. 


C 

n 


(to  /oj)‘  (u)  /w)2(v/w) 

— E - i  —2 - 

l+(v/u)2  l+(  v/uj)  2 


where  to  *  2tt t  is  the  radian  signal  frequency,  v  is  the  collision  frequency  for 
momentum  transfer,  and  the  plasma  frequency  is  given  by 

2 

op  ■  n  /m  € 
p  e  e  o 


2^Heald,  M.A.  and  C.B.  Wharton,  Plasma  Diagnostics  with  Microwaves,  John  Wiley 
&  Sons,  Inc.,  New  York,  NY,  p.  6,  1965. 
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The  permittivity  of  dree  space  is  designated  by  e  ,  the  electron  charge  by 

e,  the  electron  number  densitv  by  n  ,  and  m  is  the  electron  mass.  A  time 

e  i 

iwt  . 

ie-'ena -jnce  ot  e  .u  ;  :e  ;n  rssumeu. 


"or  a  .specified  signal  fn 

tron  density  n  tor  which  u)  *  u. 

J  ec  p 

n  <  n  and  overuense  when  n  > 
e  ec  e  — 

constant  of  unity  so  that  all 


.  1’he  plasma  is  said  to 

n  .  An  overdense  medium 
ec 

signals  ire  reflected  and 


be  underuense  vnen 
has  a  reflective 
none  transmitted 


through  the  plasma. 


Initial  estimates  of  the  effect  of  the  rocker,  plume  on  RF  propagation 

26 

from  the  missile  to  a  ground  station  were  based  on  ray  optics.  That  is, 

the  incident  signal  va  i  considered  to  be  a  bundle  of  rays  which  could  be 

individually  traced  through  the  plume.  This  type  of  analysis  is  based  on  the 

.  .  27  . 

so-called  Eikonal  equation  which  describes  the  ray  trajectories.  Absorp¬ 
tion  of  the  signal  along  each  ray  ptrh  can  be  -imputed;  by  comparing  the  diver¬ 
gence  of  the  incident  bundle  of  rays  with  that  of  the  existing  bundle  of  rays, 

the  effects  of  refraction  can  be  determined.  The  result  is  that  for  a  given 
effective  radiated  power  incident  on  Che  plasma  at  some  angle,  an  effective 
radiated  power  (ERP)  can  be  computed  at  a  somewhat  different  exiting  angle. 
Once  the  transformation  from  incident  to  exit  angle  has  been  determined,  the 
attenuation  of  the  signal  received  by  a  ground  station  can  be  obtained  by  com¬ 
paring  the  ERP  with  and  without  the  plume  present.  It  is  possible  for  an  exit 
angle  to  correspond  to  more  than  one  incidence  angle,  in  which  case  the 
signals  muot  be  summed  with  the  appropriate  phasing. 


^Malraua,  pt>  "Raybend — A  Tracing  Program  for  Microwaves  Propagating  through 
Rocket  Plumes,"  paper  presented  at  the  JANNAF  12th  Plume  Technology  Meeting, 
October  1980. 

^Born  &  Wolf,  Principles  of  Optics,  Pergamon  Press,  London,  p.  110,  1959. 


18 


It  has  been  recognized  that  the  ray  optics  treatment  tends  to  over¬ 
predict  the  attenuation  produced  fcy  the  plasma  plume,  especially  in  the  case  of 
large  attenuations.''’""  This  is  because  the  analysis  is  based  on  the  assump¬ 
tion  of  an  infinitesimal  wavelength  ind  ignores  the  signal  diffracted  around 
the  plume.  Early  investigators  of  the  RF-plume  Interaction  assumed  that  the 

diffraction  estimates  should  be  based  on  a  blockage  of  the  incident  wave  by 

’’b  1 7 

the  overdense  core  of  the  plume."  ’  ”  More  recently,  it  has  been  noted  that  a 
somewhat  larger  area  is  blocked,  since  both  the  rays  incident  from  the  trans¬ 


mitter  and  rays  projected  from  a  ground  station  are  refracted  away  from  the 

■no 

overdense  core  of  the'  plume. An  effective  blocked  area  can  be  determined  by 

using  a  ray  tracing  program  and  determining  the  intersection  of  the  rays 

approaching  closest  to  the  overdense  core.  In  these  two  diffraction  estimates, 
the  concentration  of  the  rays  at  the  edge  by  refraction  and  the  attenuation 
along  the  ray  path  are  ignored. 


An  alternate  technique  has  been  developed  by  Boynton,  Davies, 
Rajasekhar,  and  Thomson. The  effect  of  the  plume  is  analyzed  by  dissecting 
the  plasma  into  a  sequence  of  absorbing  phase  screens.  The  incident  field  is 
evaluated  on  a  planar  surface;  the  plasma  between  that  plane  and  the  succeeding 
plane  is  assumed  to  modify  the  field  at  that  plane  in  phase  and  amplitude. 


Propagation  to  the  next  plane  is  calculated  using  the  Fresnel-KLirchof  f 
integral.  Succeeding  steps  are  handled  in  the  same  way  until  the  electron 
density  has  dropped  to  a  point  that  phase  and  amplitude  changes  are  negligible. 
Fast  Fourier  transform  techniques  are  used  to  reduce  the  computation  time. 
From  the  Fourier  transform  of  the  field  at  the  last  plane,  the  far-field  radia¬ 
tion  may  be  directly  evaluated.  This  technique  includes,  in  principle,  both 
refraction  and  diffraction  and  has  the  additional  capability  that  the  gross 
features  of  plasma  turbulence  can  be  included. 


28vicente,  F.A.,  E.C.  Taylor,  and  R.W.  Phelps,  "Analysis  of  Flame  Effects  in 
Measured  Electromagnetic  Propagation  Data,"  J.  Spacecraft,  4,  pp.  1069-1075, 
1967. 

“^Peng,  S.Y.,  "Analytical  Model  for  Predicting  R-F  Amplitude  and  Phase  Effects 
of  Plume,"  JANNAF  10th  Plume  Technology  Meeting,  1977. 

J0Boynton,  F.W.,  A.R.  Davies,  P.S.  Rajasekhar,  and  J.A.  Thompson,  "Effects  of 
Large  Absorbing  Rocket  Exhaust  Plumes  on  R.F.  Signals,"  JANNAF  10th  Plume 
Technology  Meeting,  CPI  Publications  291,  pp.  143-173,  1977. 
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In  this  study,  we  shall  first  compare  the 
different  analyses — ray  optics,  diffraction,  and 
— with  the  computed  result  for  a  scattering 
u'.actly.  The  mcd  will  consist  1  radially  i 


radiated  fields  from  the  three 
the  phase  screen  approximation 
problem  which  can  be  solved 


excited  by  either  an  electric  or  magnetic  Line  source.  The  use  of  both  elec¬ 
tric  and  magnetic  line  sources  will  allow  the  assessment  of  the  effects  of  the 
two  different  polarizations  on  the  result.  The  plasma  parameters  will  be  simi¬ 
lar  to  those  computed  for  the  cnird  stage  of  the  MX  missile.  Although  the 
problem  analyzed  will  be  only  two-dimensional,  it  will  provide  a  direct  com¬ 
parison  with  the  results  of  the  diffraction  analysts  of  che  MX  plume,  since 
that  model  is  also  two-dimensional.  Comparison  will  be  made  with  the  diffrac¬ 
tion  model,  selecting  the  strip  size  based  on  the  assumption  of  an  overdense 
core  or  using  the  ray  optics  method  previously  described.  Alternate  selection 
procedures  will  aLso  be  considered. 


Similarly,  the  ray  optics  prediction  will  be  compared  with  the  exact 
solution  to  see  if  there  exists  any  region  of  validity  for  its  use,  particular¬ 
ly  at  grazing  angles  to  the  plume. 


Finally,  a  two-dimensional  formulation  of  the  absorbing  phase  screen 
model  will  be  compared  with  the  exact  solution.  Of  particular  interest  will  be 
the  manner  in  which  the  step  sizes  must  be  selected  for  a  given  signal  fre¬ 
quency  and  plasma  parameters. 

It  is  anticipated  that  the  absorbing  phase  screen  model  will  provide  the 
best  approach  to  modeling  the  interaction  of  the  RF  with  the  plasma  plume; 
however,  validation  of  this  model  is  currently  unavailable. 


After  satisfactory  comparison  is  obtained  between  the  absorbing  phase 
screen  model  and  the  exact  solution  for  a  cylindrical  plasma,  an  RF-plasma 
interaction  code  will  be  written  to  use  the  plasma  properties  generated  by  the 
plume  flow  field  codes  in  order  to  compute  the  effects  of  the  plasma  on  the 
phase  and  amplitude  of  the  RF  wave  propagating  through  it. 
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3.  RESULTS 


til  w  i  \  L  J  i  u  U  U  y  )  'i  in  i'll  13  l  j  ‘.v  u  o  p  *.  d  C  u  OH  A  -1 .11  i  »'  *  il *-  1  i  ci  iti’.U  l 

plane  density  d  i  $  t  r  ibut  i  .o  (or,  equivalently  ,  the  electron  den.; 

tion)  to  various  factors  in  the  plume  flow  field.  The  aim  is 
theoretical  plume  flow  field  model  in  which  all  the  pertinent 
included.  The  results  of  these  studies  are  given  below. 


w  i  7  L  Z.  7  01!  u  A  ^ 

ity  distribu¬ 
te  define  a 
phenomena  are 


3  .  1.  EFFECTS  OF  IONIC  REACTION  RATE 

As  mentioned  previously,  finite  rate  chemical  reactions  must  be  included 

in  the  rocket  nozzle  flow  calculation  in  order  to  determine  the  electron  mole 
fraction  at  the  nozzle  exit.  Doth  neutral  and  ionic  reactions  need  to  be  con¬ 
sidered.  The  neutral  reactions  can  affect  the  flow  field  temperature  and  the 
concentrations  of  the  neutral  species  that  enter  into  ionic  reactions.  Both 
have  effects  on  the  ionic  reaction  rate  and,  therefore,  the  electron  concentra¬ 
tion.  For  a  reaction  in  which  the  reactants  A  and  B  are  converted  into  the 
products  C  and  D 


A  +  B  *  C  +  D 
k 


the  specification  of  the  forward  reaction  rate  coefficient  k^  is  sufficient, 
since  the  reverse  rate  coefficient  k^  can  be  computed  from  the  equilibrium 
constant  K.  Thus 


k 

r 


kf/K 


(1) 


where  k^  varies  with  the  gas  temperature  T  as  follows: 

kf  -  AT_nexp(-E/RT)  (2) 

where  R  is  the  gas  constant  and  the  constants  n  and  E  are  determined  from 
experimental  reaction  rate  data. 
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For  a  propellant  combination  containing  the  chemical  elements  C,  N,  0, 

H,  Cl,  M,  K,  and  Na,  the  neutral  and  ionic  reactions  necessary  to  describe  the 

chemical  system  are  listed  in  the  Appendix  together  with  the  constants  in 

3  L 

Zq.  (2)  and  their  uncertainty  limits. 

An  examination  of  the  reaction  rates  indicates  that  large  uncertainties 
exist  in  their  values,  especially  the  ionic  reactions  which  have  the  dominant 
effect  on  the  electron  concentration.  In  order  to  evaluate  the  uncertainties 
of  the  electron  concentration  in  the  plume  due  to  the  inaccuracies  of  the  reac¬ 
tion  rates,  the  eLectron  mole  fraction  at  the  nozzle  exit  plane  was  calculated 
for  the  nominal  value  and  the  upper  and  lower  bounds  of  the  ionic  reaction 
rates  as  listed  in  Table  1  and  the  Appendix.  The  results  are  listed  in 
Table  2  for  a  solid  propellant  containing  80  ppm  sodium  and  30  ppm  potassium. 

Table  1.  Ionic  Reactions 


R19 

K  +  e  *  M  »  K  * 

M 

R20 

Na+  +  e  +  M  *  Na 

+  M 

R21 

Cl  +  e  +  M  -  Cl~ 

+  M 

R22 

Kf  Cl"  3  K  +  Cl 

R23 

Na+  +  Ci”  =  Na  + 

CJl 

R24 

HC l  +  e  *  H  +  CJL' 

- 

^ ^Jensen,  D.E.  and  G.A.  Jones,  "Gas-Phase  Reaction  Rate  Coefficients  for  Rocke¬ 
try  Applications,"  Rocket  Propulsion  Establishment  Report  No.  71/9,  Ministry 
of  Defense,  London  W.C.  2,  October  1971. 
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Table  2.  Effects  of  Ionic  Reaction  Rates  on  Electron 
Mole  Fraction  at  Nozzle  Exit 


Nominal  t 

3.0  : 

:<  10“ ; 

5  x  R19  and  R20  I 

6.27 

x  10“7 

1/5  x  R19  and  R20  j 

1.13 

x  10“° 

JO  x  R21  1 

9.2  : 

x  10"‘ 

1/30  x  R21  i 

7.43 

x  10"7 

10  x  R22  and  R23  1 

8 .46 

x  10"7 

1/100  x  R22  and  R23 

o .  88 

x  10“7 

JO  x  R24 

2.8 

✓s 

O 

1 

a 

1/100  x  R24  | 

3.86 

* 

o 

i 

It  appears  that  no  significant  change  in  the  electron  density  is  caused 

by  the  uncertainties  in  the  first  five  reactions.  The  two-body  attachment 

reaction  (.Reaction  24),  however,  has  the  predominant  effect  on  the  ionization 

level  as  shown  in  Table  2  and  Fig.  1  where  X  denotes  the  electron  mole 

e 

fraction,  R*  the  throat  radius,  and  Z  the  axial  distance  from  the  throat.  The 
uncertainties  in  its  reaction  rate  can  cause  a  factor  of  three  increase  or  a 
factor  of  two  decrease  in  the  electron  density  from  the  value  predicted  on  the 
basis  of  its  nominal  rate.  Therefore,  based  upon  the  state  of  the  art  in  the 
reaction  rates,  an  uncertainty  factor  of  about  seven  in  the  plume  electron 
density  is  to  be  expected. 

This  situation  will  prevail  until  such  time  that  the  rate  of  Reaction 
24  can  be  more  accurately  determined  by  laboratory  experiment.  In  the  absence 
of  more  accurate  reaction  rate  data,  the  recommended  approach  is  to  vary  the 
rate  of  Reaction  24  to  match  existing  flight  data,  such  as  some  of  the  S-band 
telemetry  data  available  for  Minuteman  Stage  III  flights  from  Vandenberg.  The 
same  reaction  rate  is  then  used  to  predict  plume  effects  on  RF  signals  for 
future  flights. 


iectron  Attachment  Rato 


The  dominant  effect  of  Reaction  24  can  be  traced  to  the  fact  that  in  the 


nozzle  the  gas  density  is  rapidly  decreasing.  Since  the  rate  of  Reaction  24  is 
p  report  lon.il  to  Che  density  square  while  those  of  Reactions  19  to  21  are  pro¬ 
portional  to  tne  density  cube,  the  flaw  expansion  in  ;t;  no.szie  causes  the 
Reactions  19  through  21  to  slow  down  much  more.  Reaction  2d  obviously  proceeds 
in  the  reverse  direction  in  the  nozzle  to  generate  more  electrons  from  the 
negative  chlorine  ion.  Since  the  concentration  of  atomic  hydrogen  is  orders  of 
magnitude  higher  than  the  other  reactants  in  Reactions  19  to  23,  an  uncertain¬ 
ty  in  the  rate  of  Reaction  24  produces  much  greater  error  in  the  electron 
density.  This  is  illustrated  in  Figs.  2  and  3  where  the  species  from  a  solid 
propellant  exhaust  are  plotted.  This  propellant  contains  a  low  percentage  of 
ammonium  perchlorate  so  that  the  amount  of  HC2.  formed  is  only  1  percent  by 
mole.  Also  AJl  +  exists  in  the  nozzle  but  was  not  Included  in  Fig.  3.  The 
major  species  (CO,  H, ,  H^O,  CO, ,  and  HCi)  are  essentially  frozen  during  the 
flow  expansion  in  the  nozzle.  However,  the  minor  neutral  species  and  the  ionic 
species  undergo  considerable  change.  This  demonstrates  the  importance  of 
accounting  for  finite  rate  chemical  reactions  in  the  nozzle  when  ionization  of 
the  rocket  engine  exhaust  is  considered.  The  ionization  reactions  are  so  fast 
that  considerable  care  must  be  exercised  in  choosing  a  sufficiently  small 
integration  step  size  to  integrate  the  species  conservation  equations.  It  was 
discovered  that  by  reducing  the  integration  step  size  progressively,  the 
solution  no  longer  changed  when  a  step  size  equal  to  or  less  than  0.2  percent 
of  the  throat  radius  was  taken.  This  is  shown  in  Fig.  4  together  with  a 
solution  using  a  step  size  equal  to  1  percent  of  the  throat  radius.  Clearly, 
large  errors  can  result  when  the  integration  step  size  is  too  large.  A  step 
size  of  Az/R*  *  0.002  was  adopted  for  all  chemical  kinetic  calculations  in 
this  study. 


3.2 


TWO-PHASE  FLOW  IN  ROCKET  NOZZLE 


Liquid  propellants  generally  do  not  form  condensed  solid  particles 
except  perhaps  from  erosion  of  the  nozzle  throat.  In  that  case,  the  amount  of 
particulates  carried  by  the  exhaust  gas  is  so  small  that  effects  on  the  flow 
field  can  be  neglected. 
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Solid  propellants  generally  contain  aluminum  as  part  of  the  fuel.  As 
the  propellant  burns,  the  aluminum  is  oxidized,  forming  liquid  aluminum  oxide 
particles  of  various  sizes  and  ocher  gaseous  aluminum  compounds.  The  liquid 
droplets  solidify  into  solid  particles  when  they  tool  to  the  melting  point  ;t 
aluminum  oxide  by  losing  heat  to  the  expanding  gas  in  tha  nozzle.  The  amount 
of  particulates  formed  can  be  a  large  fraction  of  the  total  exhaust,  generally 
from  20  to  40  percent,  depending  upon  the  propellant  formulation. 


The  presence  of  large  amounts  of  particulates  in  the  rocker,  nozzle 
causes  considerable  momentum  and  heat  transfer  between  the  gas  and  the 

particles.  A  fully  coupled  two-phase  flow  solution  was  obtained  for  the  MX 
third  stage  nozzle  geometry,  using  the  Aerospace  CY3D2P  code.  Perfect  gas  and 

single  size  particles  of  6.7  Um  diameter  were  assumed.  Table  3  lists  the 
gas  and  solid  properties  used  in  the  calculation. 


Table  3.,  Gas  and  Solid  Properties  for  Two-Phase 
Nozzle  Flow  Calculation 


0.1 -t  Phase 

Particle  Phase 

Ratio  of  specific  heats  «  1.31ft 

Particle  diameter  ■  6.7  pm 

Heat  capacity  «  0.436  Btu/lb-°R 

Particle  heat  capacity  •  0.326  Btu/lb-°R 

Brand ti  number  «  0,5 

Solid  mass  densltv  »  200  Ib/ft"* 

EmlssiVlty  »  0.0 

F.mlsslvlty  ■  0.25 

Chamber  pressure  ■  693  psla 

Percentage  by  weight  of  particles  ■  30. 5Z 

Chamber  tempera  tre  ■  6000°R 

Viscosity  at  6000 *R  -  3.67  x  10-5lb/ft-sec 

Temperature  exponent  for  viscosity  ■  0.65 

29 


Figure  5  shows  the  MX  third  stage  motor  configuration  and  the  domain 
used  in  the  transonic  analysis.  The  expanded  view  of  this  domain  and  the  grid 
system  tor  the  computation  are  shown  in  Fig.  o.  At  the  left  side  of  the  tran¬ 
sonic  domain,  the  gas  is  subsonic.  It  attains  sonic  velocity  in  the  throat 
region  and  passes  through  the  right-hand  computational  boundary  s’.-:ers,'-nicaLiy . 
In  the  transonic  domain,  the  unsteady  two-phase  governing  equations  are  solved 
by  finite  differences  to  avoid  the  difficulty  of  having  to  solve  mixed 
elliptical-hyperbolic  differential  equations.  The  computation  marches  in  time 
until  the  flow  variables  no  longer  change  La  time  and  the  solution  roaches 
steady  state.  The  solution  is  continued  into  the  expanding  nozzle  until  the 
nozzle  exit  is  reached.  In  this  region,  the  steady-state  governing  equations 
are  solved.  The  computational  grids  at  a  number  of  selected  axial  stations 
are  depicted  in  Fig.  7. 


The  results  of  the  comput.it ion  with  and  without  the  particles  are  shewn 
in  Figs.  8  through  15  to  demonstrate  the  effects  of  the  particles  on  the  nozzle 
flow.  Contours  of  constant  Mach  number  and  ratio  of  gas  density  to  chamber 
density  from  the  two-phase  flow  solution  aie  plotted  in  Figs.  8  through  11. 
The  large  radial  gradient  of  the  Mach  number  near  the  nozzle  wall  as  shown  in 
Fig.  8  is  the  result  of  the  absence  of  the  particulates.  The  particles  cannot 
negotiate  the  large  throat  curvature,  and  a  limiting  particle  streamline  is 
formed.  At  nozzle  radius  larger  than  the  boundary  of  this  limiting  streamline, 
the  region  is  occupied  only  by  the  gas. 

Similar  contour  plots  from  the  solution  without  the  particles  are  pre¬ 
sented  in  Figs.  12  through  15.  Comparison  of  the  results  with  and  without  the 
particles  shows  that  the  presence  of  the  particles  greatly  decreases  the  gas 
Mach  number,  while  the  gas  density  is  relatively  unaffected. 


In  order  to  check  the  validity  of  the  SPP  code,  the  same  case  was  com¬ 
puted  using  identical  gas  and  particle  properties  and  the  perfect  gas  mode  of 
the  SPP  code.  Radial  profiles  of  gas  properties  at  the  nozzle  exit  plane  are 
compared  with  those  obtained  from  the  Aerospace  CY3D2P  code  in  Figs.  16  and  17. 
Satisfactory  agreement  of  the  results  from  the  two  codes  is  obtained  for  Mach 
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Fig.  5.  MX  Third  Stage  Motor  Interior  Configuration 


Fig.  6.  Computational  Grid  for  Transonic  Flow  Regime 
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Fig.  8.  Mac'n  Number  Contour,  Two-Fhaae 
(Transonic  Regime) 


Gas  Density  Contour,  Two-Phase 
(Transonic  Regime) 


Fig.  10.  Mach  Number  Contour,  Two-Phase 
(Supersonic  Regime) 


Fig.  11.  Gas  Density  Contour,  Two-Phase 
(Supersonic  Regime) 
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Fig.  12.  Mach  Number  Contour,  One-Phase 
(Transonic  Regime) 


Fig.  13.  Gas  Density  Contour,  One-Phase 
(Transonic  Regime) 
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Fig.  14.  Mach  Number  Contour,  One-Phas 
(Supersonic  Regime) 


Fig.  15.  Gas  Density  Contour,  One-Phase 
(Supersonic  Regime) 
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Profiles  of  Pressure  and  Density  at  Nozzle  ExiL  Plane 
with  Single  Particle  Size  of  6.7  |im 


number,  temperature,  pressure,  and  density.  The  small  differences  can  be 
explained  in  part  by  the  different  particle  drag  and  heat  transfer  coefficients 

adopted  by  the  two  codes.  It  is  interesting  to  note  that  the  >as  temperature 

13  auwUC  ’  a  II  1  (1  ti  1“  W  ft  d  f\  pUlTCioittS  3  £'3  p  T  3  \~l  C  »  l  m  p  ii  T  £  1 C  i  3  ,  0  J  O  ..  1 1 1  £  -t  C  .1 

lower  rate,  neat  up  the  gas.  As  a  result,  the  gas  temperature  is  increased  in 
the  core.  The  particles  are  absent  near  the  nozzle  wall,  and  the  gas  tempera¬ 
ture  starts  to  approach  the  value  for  one-dimensional  gas  flow  without  parti¬ 
cles.  The  lowering  of  the  Mach  number  in  the  two-phase  flow  is  largely  due  to 
the  higher  gas  temperature  as  the  sonic  velocity  is  greater. 

With  the  SPP  code  validated  for  the  case  of  a  single-size  particle,  two- 
phase  flow  computation  for  a  distribution  of  particle  sizes  was  made.  For  a 
4.32-in.  throat  radius  of  the  MX  third  stage,  the  SPP  code  predicts  a  mean 
particle  diameter  of  6.7  lim.  A  log  normal  distribution  is  then  assume!  to 
divide  the  particles  into  five  diameters  of  equal  weight,  viz.,  3.01,  4.89, 
6.34,  9.36,  and  15.6  dm.  Figure  13  shows  the  gas  temperature  and  density 
profiles  at  the  nozzle  .exit  plane.  The  gas  temperature  increase  in  the  core 
is  again  seen  when  particles  are  present.  The  radial  temperature  gradient  is 
much  steeper  in  Fig.  18  than  that  shown  in  Fig.  16.  This  is  due  to  the  multi¬ 
ple  particle  sizes.  The  drag-to-weight  ratio  of  the  particles  is  smaller  as 
the  particle  diameter  is  increased.  Consequently,  the  gas  cannot  turn  the 
larger  particles  sufficiently  to  reach  the  nozzle  lip.  The  larger  particles, 
therefore,  remain  in  the  core.  The  limiting  streamline  locations  shown  in 
Fig.  18  for  the  various  particle  sizes  define  the  largest  radial  location  in 
the  nozzle  such  particles  can  extend.  Thus,  more  particles  are  located  closer 
to  the  nozzle  axis  and  cause  more  heating  of  the  gas.  The  shape  of  the  gas 
temperature  profile  is  the  result  of  this  separation  of  particle  sizes. 

These  results  for  the  case  of  multiple  particle  sizes  differ  from  those 
of  single-size  particles  not  only  because  of  the  particle  size  effect  but  also 
because  of  the  differences  in  the  chamber  condition  and  the  gas  properties  in 
the  two  computations. 


39 


ARROWS  SHOW  LIMITING  STREAMLINE 


Figure  19  shows  Che  Mach  number  M  and  pressure  p  variation  across  Che 
nozzle  exit.  The  lowering  of  the  Mach  number  in  Che  core  is  undoubtedly 
related  to  the  higher  gas  temperatures.  However,  the  gas  also  loses  consid¬ 
erable  momentum  in  acceieratins  the  particulates  as  shown  by  the  product 
( 1/2 )'  pM~ ,  the  gas  momentum  flux. 

The  mechanisms  for  creating  particles  of  various  sizes  are  complex  and 
not  well  understood.  The  combustion  process  undoubtedly  plays  a  part.  It  was 
discovered,  however,  that  the  mean  particle  size  is  somewhat  dependent  upon 
the  threat  radius.  As  the  particles  travel  through  the  throat  region,  they 
are  still  in  a  molten  form.  The  shear  force  in  the  throat  region  causes  the 
particles  to  break  up  into  smaller  sizes.  Some  correlation  of  experimental 

data  does  show  the  mean  particle  size  to  increase  with  the  throat  radius  (less 

3? 

shear  forces).  A.  correlation  by  Radke,  et  al.  "  yields  9  Mm  rae.nr  particle 
diameter  for  the  nozzle,  throat  radius  of  the  MX  third  stage  engine.  In  order 
to  determine  the  sensitivity  of  the  particle  sizes,  the  two-phase  nozzle  flow 
analysis  was  repeated  for  the  increase  in  the  mean  particle  diameter  from 
6.7  Mm  to  9  Mm.  Comparison  of  the  exit  Mach  number  profiles  in  Fig.  20 
shows  that  the  uncertainty  in  the  particle  size  does  not  introduce  large 
deviations  in  the  gas  flow  properties,  unless  the  mean  particle  diameter  is 
grossly  in  error. 

We  have  shown  that  the  presence  of  the  particulates  in  the  nozzle 
causes  appreciable  changes  in  the  gas  properties  at  the  nozzle  exit  plane. 
This  wili  subsequently  affect  the  gas  density  (and  electron  density)  in  the 
exhaust  plume,  the  extent  of  which  will  be  discussed  in  Section  3.3.2. 


32padke,  H.H.,  L.J.  Delaney,  and  P.  oraith,  "Exhaust  Particle  Size  Data  from 
Small  and  Large  Solid  Rocket  Motors,"  Report  No.  TOR-1001(S2951-l8)-3 ,  The 
Aerospace  Corporation,  El  Segundo,  CA,  July  1967. 
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3.3 


ROCKET  EXHAUST  PLUME  FLOW  FIELD 


With  the  i  Low  properties  determined  ac  che  no:zle  exit  p  1  me ,  the  c  jmpu- 
tit  ion  can  Chen  into  eke  jxhaust  plume.  At  very  high  altitudes,  v.v.ar.s 
the  atmosphere  ia  vary  tain,  the  rocket  exhaust  can  be  assumed  to  expand  into 
a  vacuum  and  che  gas  density  decrease  raonotonical ly  with  increasing  discance 
from  the  nozzle  and  with  increasing  polar  angle  from  the  plume  centerline.  At 
lower  altitudes,  however,  interaction  between  the  atmosphere  and  the  plume  g  is 
confines  Che  plume  into  a  region  the  boundary  of  which  is  determined  by  a 
balance  of  che  dynamic  pressures  of  the  two  mixing  gases. 


A  study  was  made  to  determine  the  effects  of  various  physical  phenomena 

on  the  plume  gas  density  and  electron  density  distributions  in  order  to  define 
a  proper  plume  model  in  which  the  significant  physical  phenomena  are  included. 
The  results  are  presented  below. 


3.3.1  Approximate  Far-Field  Plume  Models 

A  number  of  serai-empirical  plume  models  are  available  to  compute  the 
plume  flow  field  when  atmospheric  interaction  can  be  neglected.  These  approxi¬ 
mate  models  have  been  successfully  applied  to  the  far  field  of  the  plume  but 
become  inaccurate  as  the  nozzle  exit  is  approached  and  the  near  field  expansion 
of  the  nozzle  flow  is  felt. 


Two  such  models  were  investigated  for  the  MX  third  stage  vacuum  plume 

17  18 

application,  viz.,  the  Brook  model  and  the  Boynton  model.  The  former 
model,  however,  requires  some  adjustment  of  parameters  to  match  an  exact  solu¬ 
tion  so  that  it  is  difficult  to  apply  from  knowledge  of  the  rocket  engine 
characteristics  alone. 


with 


In  the  Brook  model,  the  plume  density  is  computed  from 

pVr^  *  F(r,m)exp  -A^(L-cos0)m 

p  =  p/p  ,  V  =  V/V*  ,  r  =  r/R 
o  e 


(3) 

(4) 
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i  whe  re 

3  .!  r.  i  ;nat  ior  lens  icy 

3  sonic  velocity  at  the  nozzle  threat 
Re  3  nozzLe  exit  radius 
r  3  radial  distance  from  the  nozzle 
0  *  polar  angle 

The  A  depends  on  the  rocket  engine  and  nozzle  characteristics,  and  m  is  a 

parameter  to  be  matched  to  an  exact  solution.  Along  the  plume  axis,  6*0 

and  Eq.  (1)  becomes 

pVr"  *  Rr,ffli  (5) 


From  a  method  of  characteristics  soLution  for  uniform  flow  properties  at 
the  nozzle  exit  plane  Oi.e.,  with  the  particulate  effects  in  the  nozzle  flow 
neglected),  the  function  F,  which  represents  the  mass  flow  per  unit  solid  angle 
along  the  plume  axis,  is  evaluated  and  the  parameter  m  in  Eq.  (5)  is  adjusted 
to  yield  the  best  fit  (see  Fig.  21).  Note  that  m  »  0.825  provides  the  best  fit 
to  the  exact  solution  in  the  far  field,  and  large  errors  are  expected  for  small 
r.  With  the  parameter  ra  fixed,  the  plume  density  can  be  computed  from  Eq.  (3) 
and  any  other  flow  properties  computed  from  isentropic  relations. 

The  Boynton  model  gives  the  following  expression  for  the  plume  density: 


with 


(7) 
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wiie  re 


3  limitin'.;  {.l.i  velocity 
n 


,  3  maximum  ?r-uvi  c  L-M*i”«r  cumins  -mule 
p*  3  density  at  the  nozzle  throat 
V*  *  sonic  velocity  at  the  nozzle  throat 
R*  3  throat  radius 

The  constant  C  is  close  to  1,  and  for  the  present  case  C  3  0.616. 

In  Fig.  22,  predictions  of  the  constant  electron  density  contours  from 

the  two  approximate  models  ate  compared  to  the  exact  solutiin  >r  th<’  '-IX  third 
stage  plume  with  uniform  nozzle  exit  properties  and  exit  electron  mole  fraction 

-7  9 

of  8.6  x  10  .  The  contours  labeled  ng  *  2.156  x  10  correspond  to  the  over- 

dense  boundary  for  the  UHF  command-de struct  frequency  of  416.5  MHz.  It  is 
evident  that  the  Brook  model  yields  better  agreement  with  the  exact  solution. 

These  studies  indicate  that  for  liquid  rocket  engines  where  the  amount 
of  particulate  in  the  exhaust  is  nil  or  insignificant,  the  approximate  vacuum 
plume  models  will  provide  quick  and  simple  estimates  of  the  plume  electron 
density  distributions.  However,  when  large  amounts  of  particulates  are  pres¬ 
ent,  as  in  Pne  caae  of  solid  rocket  engines,  the  nozzle  exit  flow  properties 
are  nonuniform.  They  cannot  be  accounted  for  by  the  approximate  models,  and 
numerical  solutions  muse  be  found. 

3.3.2  Effect  of  Particulates  on  Plume  Electron  Density 


A  study  was  made  to  evaluate  the  effect  of  the  presence  of  solid  parti¬ 
cles  in  the  MX  third  stage  plume  on  the  gas  density  and  electron  density  dis¬ 
tributions.  The  nozzle  exit  plane  flow  properties  as  shown  in  Figs.  18  and  19 
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Contours  of  Constant  Electron  Density  for  Uniform  lixiL  Mach  Number 


were  used  to  continue  the  plume  flow  field  calculation  utilizing  both  the  TD2P 
and  the  MOC  code.  The  latter  code  solves  the  inviscid  equations  of  motion  for 


o  n  l  v . 


:  cover tea 


d  e  n  s  i  t  v  c  o  n  1 3  u  r  s 


Apparently,  the  presence  of  30  percent  A particles  in  the 
plume  ha3  only  a  small  effect  on  the  electron  density  distribution.  The  par¬ 
ticulates  cause  large  changes  in  the  gas  properties  of  the  rocket  nozzle,  as 
shown  previously,  because  of  exchange  of  heat  and  momentum  between  the  two 
phases.  This  effect  was  included  in  all  the  results  presented  in  Fig.  23. 
When  the  rocket  exhaust  expands  into  a  much  lower  pressure  outside  the  nozzle, 
tne  gas  density  is  greatly  reduced.  As  a  result,  the  heat  and  momentum  trans¬ 
fer  between  the  two  phases  quickly  ceases,  and  the  gas  and  particle  motions 
become  essentially  independent  of  each  other.  The  close  agreement  between  the 
gas-only  and  the  two-phase  flow  solutions  demonstrated  this  decoupling  of  the 
motions  of  the  two  phases. 


As  the  two-phase  plume  calculation  is  more  complex  and  time-consuming 
than  the  gas-only  calculation,  the  latter  approach  appears  more  attractive  and 
economical  and  should  generate  electron  density  distribution  of  sufficient 
accuracy  for  the  subsequent  RF  interaction  calculation.  For  these  reasons, 
results  shown  in  the  remainder  of  this  report  were  generated  for  the  gas  phase 
plume  only. 


However,  the  presence  of  the  particulates  affects  the  plume  density  and 
electron  density  distributions  through  their  interaction  with  the  gas  phase  in 
the  rocket  nozzle  to  produce  nonuniform  nozzle  exit  properties  as  shown  previ¬ 
ously.  Thi3  is  shown  in  Fig.  24  where  gas  plume  calculations  were  made  at  the 
310  kft  altitude,  using  both  the  constant  and  the  nonuniforra  flow  properties 
in  Figs.  18  and  19  as  inputs.  The  latter  properties  reflect  the  effects  of  the 
particulates.  At  the  310  kft  altitude,  the  atmosphere  is  sufficiently  dense 
to  confine  the  plume  to  within  the  plume  boundary  as  shown  in  Fig.  24  with  the 
generation  of  the  barrel  shock  interior  to  the  boundary.  Figure  24  shows  that 
the  presence  of  the  particulates  in  the  nozzle  increases  the  plume  size  and 
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Electron  Density  Contours  for  Gas  and  Two— Phase  Flume 


PLUME  BOUNDAR 


causes  significant  changes  in  the  gas  density  and  electron  density  distribu¬ 
tions.  For  solid  rocket  engines,  therefore,  a  two-phase  flow  calculation  in 
the  rocket  nozzle  should  be  made  to  account  for  the  effects  of  the 
part icul -ices . 

J .  .3 . 3  effect  of  external  Atmosphere 

When  the  ignition  altitude  of  the  rocket  engine  is  sufficiently  low, 
such  as  for  the  BOA-3  trajectory  of  the  planned  MX  test  flight  where  the  third 
stage  ignites  at  259  kr't,  interaction  with  the  surrounding  atmosphere  causes 
the  plume  to  be  confined.  The  plume  boundary  will  increase  with  altitude  until 

the  atmospheric  density  is  so  small  that  the  plume  can  be  considered  to  be 

expanding  into  a  vacuum.  Figure  25  demonstrates  this  atmospheric  effect.  The 
shaded  portion  depicts  the  overdense  region  at  the  UHF  command-des truct  fre¬ 
quency  where  RF  signals  are  totally  reflected.  In  addition  to  the  overdense 
core,  there  exists  an  overdense  region  between  the  plume  boundary  and  the 
barrel  shock  due  to  the  gas  density  increase  across  the  barrel  shock.  The 
maximum  diameter  of  the  shock-layer  overdense  region  increases  with  altitude 
and  then  decreases  until  at  some  altitude  the  shock  layer  is  no  longer  over- 
dense  and  only  the  overdense  core  remains.  For  the  MX  third  stage  this  is 
estimated  to  be  at  about  310  kft.  Below  this  altitude,  the  atmospheric  effect 
should  be  considered.  However,  this  effect  occurs  only  within  a  small  portion 
of  the  entire  ignition  altitude  region.  Its  impact  on  the  overall  range  safety 
operation  needs  more  careful  evaluation. 

3.3.4  Effect  of  Nozzle  Wall  Boundary  Layer 

As  the  rocket  exhaust  flows  through  the  nozzle,  a  thin  boundary  layer 
near  the  nozzle  wall  is  generated.  For  the  large  thrust  levels  generally 
employed  for  the  various  stages  of  the  launch  vehicle,  the  boundary  layer  is 
usually  turbulent.  The  Mach  number  in  the  boundary  layer  decreases  from  super¬ 
sonic  at  the  boundary  lay  r  edge  to  subsonic  near  the  nozzle  wall.  Because  of 
the  lower  Mach  number  and  its  proximity  to  the  wall,  the  gas  in  the  boundary 
layer  expands  to  much  .larger  angles  than  the  rest  of  the  exhaust  gas  once  it 
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exits  from  the  nozzle.  The  presence  of  the  boundary  layer,  however,  is  not 
expected  to  alter  the  electron  density  significantly  in  those  regions  of  the 
plume  where  37  interne ‘.ion  is  appreciable,  since  the  electron  density  in  that 
-M  rt  of  the  plume  which  originates  from  the  boundary  layer  is  much  lower  t'vn 
in  the  interior  of  the  plume  because  of  the  large  angle  expansion  of  the  bound¬ 
ary  layer  gas.  Nevertheless,  the  effect  of  the  boundary  layer  was  investigated 
by  computing  the  inviscid  expansion  of  the  boundary  layer  with  the  MOC  code. 

In  Fig.  26,  results  with  and  without  the  presence  of  the  nozzle  boundary 
Layer  are  compared.  It  is  evident  that  the  boundary  layer  effect  is  insignifi¬ 
cant.  Although  not  shown  in  the  figure,  the  electron  density  distributions  for 

the  two  cases  are  nearly  identical.  It  is  concluded  that  for  RF-plume  inter¬ 
action  analysis  the  presence  of  the  nozzle  boundary  layer  can  be  neglected. 


3.3.5  Effect  of  Nozzle  Length 


In  the  study  of  the  MX  third  stage  pLurne  effects  on  RF  signal  propaga¬ 
tion,  the  effect  of  nozzle  extension  was  considered.  In  the  stowed  configura¬ 
tion  before  the  second  and  third  stage  separation,  the  third  stage  nozzle 
measures  36.5  in.  from  the  throat.  The  nozzle  is  extended  to  75.0  in.  when 
the  third  stage  is  deployed.  The  plume  flow  field  was  computed  for  both 
nozzle  lengths  to  determine  the  difference  in  the  election  density  distribu¬ 
tion.  The  results  are  presented  in  Fig.  27.  The  exhaust  plume  from  the  un¬ 
extended  nozzle  has  a  smaller  UHF  overdense  core.  This  is  mainly  due  to  the 
lower  electron  mole  fraction  in  the  exhaust  gas  from  the  unextended  nozzle  as 
shown  by  the  electron  mole  fraction  variation  along  the  nozzle  in  Fig.  4. 
Some  displacements  of  the  plume  boundary  and  the  barrel  shock  are  also  noted. 
However,  the  change  in  the  plume  shape  is  small. 


These  results  imply  that  the  exhaust  plume  from  the  fully  extended 
nozzle  will  cause  greater  RF  inceraction  and  should  be  chosen  as  the  more 
conservative  case. 
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=  4.538  WITH  BOUNDARY  LAYL:R 


Fig.  26.  Effect  of  Nozzle  Boundary  Layer 


4.  CONCLUSIONS 


p  ’ :  r  —  o  3  '  or  1 a £  a  ^ «. >•»  £  w  ^  <1 J  ^ a ^ *  5 


it  axhau; 


o Lu;no  or. 


RF  signal  propagation,  it  is  necessary  co  provide  plume  flow  field  properties, 
such  as  gas  density,  electron  density,  and  temperature  distributions,  from 
theoretical  prediction.  From  studies  of  the  various  factors  which  influence 
the  ionization  and  gas  density  distributions  in  the  exhaust  plume  of  a  solid 
rocket  engine,  the  following  conclusions  are  obtained: 


o  Electrons  in  the  rocket  exhaust  are  generated  by  ionization  of 
the  alkali  contaminants  in  the  rocket  chamber.  The  ionization 

level  changes  as  the  exhaust  flows  through  the  nozzle  due  to 
additional  reactions  among  the  ionic  species,  with  the  two-body 
attachment  reaction  HCfi,  +  e  ■  H  +  C2,  being  predominant  in 
controlling  the  ionization  level.  Based  upon  the  present  know¬ 
ledge  of  its  reaction  rate,  a  factor  of  seven  uncertainty  in  the 
plume  electron  density  distribution  can  be  expected. 

o  The  particulate  matters  in  the  two-phase  exhaust  affect  the  plume 
electron  density  distribution  mainly  through  their  effects  on  the 
nozzle  flow  field  as  they  interact  with  the  gas  phase  in  the 
exchange  of  heat  and  momentum.  Their  presence  in  the  plume  flow 
field  can  be  neglected  in  the  determination  of  electron  density 
distribution  due  to  the  weak  interaction  between  the  particulates 
and  the  plume  gas. 

o  In  the  altitude  range  below  310  kft,  the  plume  shock  formed  by 

the  interaction  of  the  MX  third  stage  exhaust  with  the  sur¬ 
roundings  causes  the  appearance  of  an  outer  UHF  overdense  region 
which  should  be  considered  in  the  calculation  of  RF  signal  degra¬ 
dation. 

o  The  presence  of  the  nozzle  wall  boundary  layer  has  a  negligible 

effect  on  the  plume  plasma  properties  and  the  RF  signal  propaga¬ 
tion. 
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o 


The  exhaust-,  from  the  fully  extended  nozzle  of  the  MX  third  stage 
presents  a  more  severe  plasma  environment  to  RF  signal,  propaga¬ 
tion  than  that  of  the  unextended  nozzle. 

o  Approximate  far  field  plume  models,  such  as  she  3 rook  model  and 

the  Boynton  model,  may  be  useful  for  the  description  of  plumes 
from  liquid  rocket  engines.  They  cannot  account  for  the  inter¬ 
action  between  the  particulates  and  the  gas  in  the  rocket  nozzle 
and  are  noc  accurate  enough  to  describe  plumes  from  solid  rocket 
engines. 
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APPENDIX 

CHEMICAL  REACTIONS  IN  ROCKETRY  APPLICATION 

L units) 


E 

!J  n  c  e  r  c  a  i 

.  n  c  v  F  a  e  c  0 

No. 

Reac cions 

A 

n 

Ckcal/ 
mole ) 

Upper 

Bound* 

Lowe  r 

Bound** 

1. 

M  +  H  +  OH  -  H20  +  M 

7.5  x  1023 

2.6 

0. 

10 

10 

2 . 

M  +  H  +  H  *  +  M 

6.4  x  10 17 

L.O 

0. 

30 

30 

3. 

M  +  CO  +  0  -  C02  +  M 

1.38  x  1034 

5.24 

5.6 

3 

30 

4. 

M  +  H  +  cl  -  UC2,  +  M 

3.0  x  10 16 

0.5 

0. 

100 

100 

5. 

M  +  MCI  +  Cl  -  AAC2,2  " 

•  M  3.0  x  10 16 

0.5 

0. 

LOO 

LOO 

6. 

M  +  MC%2  +  Cfi.  *  AACS-3 

+  M  3.0  x  1016 

0.5 

0. 

100 

LOO 

7. 

Ho  +  OH  -  H  +  H20  _ 

2.19  x  10u 

0. 

5.  15 

3 

3 

8. 

CO  +  OH  -  C02  +  H 

5.6  xlO11 

0. 

1.08 

5 

5 

9. 

HCJl  +  OH  *  H20  +  Ql 

10U 

-0.5 

6.0 

30 

30 

10. 

AACA2  +  H  -  AJICJl  +  HCA 

1.05  x  10u 

-0.5 

2.566 

11. 

Cfi.2  +  H  -  HC2-  +  Cl 

u 

0 

X 

0 

0. 

3.0 

5 

5 

12. 

CO  +  0  ■  C02 

1.78  x  1010 

0. 

2.53 

13. 

M  +  HCS.  -  A2.C2.  +  H 

5.0  x  1011 

-0.5 

5.673 

14. 

AJlO  +  HC2,  *  A5.0CS.  +  H 

10u 

-0.5 

5.673 

15. 

MCI  +  OH  -  AJLOCi  +  H 

10u 

-0.5 

5.619 

16. 

AJIOH  +  OH  *  H  +  HAS,02 

10U 

-0.5 

5.627 

17. 

K  +  HCJl  3  KCJl  +  H 

3.612  x  1014 

0. 

5.0 

10 

30 

18. 

Na  +  HC2.  -  NaCS.  +  H 

3.01  x  IQ14 

0. 

8.0 

10 

30 

19. 

M+K++e-K+M 

7.25  x  1025 

1.5 

0. 

5 

5 

* 

** 

Factor  by  which  the  rate 
Factor  by  which  the  rate 

is  increased 
is  decreased 
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Uncertainty  Factor 


E 


No . 

Reac  t ions 

A 

n 

(kcal/ 
mo  1  e  ) 

.  Upper 
Bound* 

Lower 

Sound' 

20. 

M  +  Ma*  +■  e  =  Na  M 

5 ,  a  3  5 

x  10  I? 

2.0 

0 . 

5 

5 

21. 

M  +  Ci.  +  a  »  Ci"  +  M 

1.09 

x  1G13 

0. 

0. 

30 

30 

->  n 

*.  « 

*  Ci"  *  ;<  +  ci 

6.02 

x  10  15 

0.3 

0. 

10 

100 

23. 

Na+  +  CZ~  =  Na  +  Ci 

1.31 

x  10  16 

0.5 

0. 

10 

100 

24 . 

HCi  f  e  -  11  +  Cl~ 

6.02 

x  10 ‘5 

0. 

19.67 

30 

100 

The  symbol  M  in  the 

three-body 

reac  t ions 

denote 

s  any  ch 

amical 

spec ies 

present  in  the  system.  The  species  have  different  catalytic  efficiencies  of 
causing  the  reaction  to  proceed  during  collisions.  The  factors  listed  below 
for  Reactions  1  and  2  are  used  to  multiply  the  individual  concentration  before 
their  concentrations  are  summed  in  the  reaction  rate  calculation.  For  all 
other  reactions,  a  catalytic  efficiency  of  1  is  used. 


Catalytic  Efficiencies  for  Three-Body  Reactions 


Spec  ies 


Reactions 

HjO 

H2 

H 

0 

OH 

CO 

£22 

l 

20 

5 

12.5 

12.5 

12.5 

8 

8 

2 

20 

4 

25.0 

25.0 

25.0 

4 

4 
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